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ABSTRACT 

An experiment was performed to confirm the 
proportionality of log-amplitude variance and the 7/6 power 
of wave number predicted by Tatarski for horizontal 
propagation from a spherical-wave transmitter to a point 
detector. The validity of this proportionality was tested for 
two wavelengths, 0.632 and 10.6 microns. A 0.632p He-Ne 
laser and a 10. 6ju C0 2 laser transmitter were operated 
simultaneously over a folded 1.1 -km horizontal path, and the 
beams were detected by a photomultiplier and a gold-doped 
germanium detector. The primary scintillation statistic, 
log-amplitude variance, was evaluated for each wavelength on 
a digital computer. The ratio of variances at 0.632 and 10.6 
microns was in close agreement with predictions. Power 
spectral density, autocorrelation, and cumulative probability 
density also were evaluated for each wavelength. Scintillation 
statistics at 10.6p were found to be log-normal, just as in the 
visible portion of the spectrum. 
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WAVELENGTH DEPENDENCE OF 
LASER BEAM SCINTILLATION* 


by 

M. W. Fitzmaurice, J. L. Bufton, and P. O. Minott 
Goddard Space Flight Center 


INTRODUCTION 

Theoretical analyses by several investigators have led to predictions that the strength of optical 
scintillation depends strongly on wavelength. Recent statistical studies are couched in terms of log 
amplitude, HO, defined as 


m = 1/2 iog 4 i(t)n 0 , (1) 

where I(t ) is the instantaneous irradiance, and 7 0 is the mean irradiance. Using the Kolmogoroff model 

of turbulence, it can be shown (References 1 and 2) that, for horizontal paths, the variance of HO (a 2 ) 
is 7 


cr* ** k 1 ! 6 z ] i/e c 2 n 


where 

k = wave number, 
z = range, 

C 2 n = refractive index structure constant. 


( 2 ) 


We have conducted an experiment to test the validity of the wavelength dependence predicted by this 
expression. 


EXPERIMENTAL TECHNIQUE 

Because of the short-term stationarity and the lack of statistical homogeneity of the atmosphere, 
the experimenter must transmit simultaneously through the same portion of the atmosphere at each 
wavelength of interest. The experiment for this report is shown schematically in Figure 1 . A He-Ne laser 
at 0.632 micron and a C0 2 laser at 10.6 microns were aligned parallel, with an offset of 12 cm; and were 
operated over a folded horizontal 1.2-km path. The He-Ne output was unmodulated and no external 
optics were used, so that the assumptions of a spherical wave transmitter were satisfied. The C0 2 output 


* Presented to the Spring 1968 Meeting of the Optical Society of America, Washington, D.C. 
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APERTURE STOP 



Figure 1— Experiment schematic diagram. 


was modulated at a 1-kHz rate to minimize the 
background and to enhance the detector operation. 
The 10.6-micron beam was transmitted through a 
13-mm-aperture telescope which, because of the 
longer wavelength, also approximated a spherical 
wave source. A small portion of the C0 2 beam was 
diverted and detected by a thermopile to monitor 
any fluctuations in output power. 

A mirror in the field reflected the two beams, 
and they were then detected at the transmitter 
station by a photomultiplier and a gold-doped 
germanium detector. The output of the 
photomultiplier tube was recorded directly on an 
FM tape recorder. The 10.6-micron detected signal 
was demodulated by means of a standard envelope 
detection scheme and then recorded. 


Table 1 lists the important parameters of the experiment. In order to avoid the effects of 
aperture-averaging, each detector aperture was made to approximate a point detector. In addition, the 
photomultiplier tube aperture was adjusted so as to make the ratio of the apertures equal to the square root 
of the ratio of the wavelengths. This ensured that aperture-averaging effects were insignificant. 


The recorded scintillation data were replayed through an A/D converter into a Raytheon 520 digital 
computer that analyzed the statistics of the scintillation. The computer was programmed to analyze both 
irradiance and log-amplitude statistics. Available statistics and functional relationships are given in Table 2. 
The analog scintillation signal was sampled at a 
1-kHz rate for 10 seconds, with each sample 
assigned to one of 500 discrete levels available in 
the computer. The selection of the time-record 
length for a single data run is a tradeoff between 
the statistical reliability of a measurement and the 
statistically time-varying atmosphere. After the 
bandwidth of the scintillation was evaluated, a 
10-second record appeared to be a near-optimum 
choice. 


RESULTS 

Figure 2 shows the results of measurements 
taken over a 2-hour nighttime period in December. 
Plotted along the Y-axis is the parameter \p defined 
as the ratio of log-amplitude variances at 0.632 and 
10.6 microns. The X-axis is the log-amplitude 
variance for the He-Ne beam and the refractive 
index structure constant C N 2 , which in this case 


Table 1 

Experimental Parameters 


Parameter 

co 2 

He-Ne 

Output power 

2 watts 

1 mw 

Beam divergence 

0.4 mr 

0.7 mr 

Beam diameter 

1 3 mm 

1 .4 mm 

Detector aperture 

2.5 cm 

.62 cm 

Detector field of view 

2 

5 

Detection signal bandwidth 

220 Hz 

2.5 kHz 

Chopper frequency 

1 kHz 


Range 

1 .2 km 

1 .2 km 

Average beam height 

4m 

4m 

Folding mirror diameter 

50.8 cm 

50.8 cm 

Length of data run 

1 0 sec 

10 sec 
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was determined from the equations for horizontal 
spherical wave propagation. The value 26.7 is the 
theoretical ratio obtained from the 7/6 power 
predictions. 


Table 2 

Computer Analysis 


The data would seem to indicate that the 
scintillation statistics depart from the 7/6 power 
law in the region to the left. The difficulty here is 
that region to the left corresponds to very low 
turbulence. The C0 2 laser output unfortunately 
had an appreciable amount of 120-Hz power 
supply ripple. In this region, the dominant noise on 
the 10.6-micron signal was caused not by the 
atmosphere but by the 120-Hz ripple. The output 
of the C0 2 unit was recorded on tape; and an 
evaluation of the magnitude of the ripple showed 
that, in the 0.04 region, the internal laser 
modulation became the same order of magnitude 
as the atmospheric-induced modulation. Therefore, 

the data points in the region below 0.07 are of questionable value. Some compensation for the C0 2 
internal noise can be accomplished if one makes appropriate assumptions. 

The intensity at the detector aperture can be considered (Reference 3) as the product of two random 
variables, 


Irradiance statistics 

Log-amplitude statistics 

Average value 

Average value 

Variance 

Variance 

Median 

Median 

3rd central moment 

3rd central moment 

4th central moment 

4th central moment 

Probability density 

Probability density 

function 

function 

Power spectral density 
Autocorrelation function 

Cumulative distribution 
function 


lit) = I L hit) 


(3) 


The first, I L , is the irradiance fluctuation caused by internal laser ripple; the second, h(f), is the 
time-varying atmospheric transmission. If both variables are assumed log-normal with respective variances 

a j 2 and a 2 2 , their product is also log-normal with 
a variance equal to a t 2 +o 2 2 . Several experiments 
(References 3 and 4) have shown hit) to be 
log-normal (in the visible spectrum), and inspection 
. _ • •. • of cumulative density of log amplitude for the 

7 .-'-v transmitter noise (Figure 3) has indicated 

._ • normality as a reasonable approximation. 


26 7 ■ 
! 25 ■ 




15h 


101 i I 
0 00 o 04 0 08 


2 8 


0.12 0 16 0 20 0.24 0.28 

cr 2 

2 £0 632 

I j ! ! 

14 20 26 32 

C N x 10 15 (meters) ~ 2 ' 3 


The data are shown replotted in Figure 4 with 
the transmitter noise subtracted from the C0 2 
statistics; the measurements obtained during 
low-scintillation conditions are neglected. The 
average value of the 39 independent data points 
shown is 26.8, and their variance is 10.8. The 
average value of i p is in close agreement with the 
predicted value of 26.7. 


Figure 2— Ratio of log-amplitude variance for He-Ne and 
C0 2 wavelengths. 


The statistical reliability of the measurements 
has been analyzed to determine the magnitude of 
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Figure 3-C0 2 transmitter noise; log amplitude-cumula- 
tive probability. 


Figure 4— Ratio of log-amplitude variances for He-Ne 

and C0 2 wavelengths (transmitter noise removed, and 

all data less than 0.07 neglected). 

j2 (0.632) 


-• 


data scatter to be expected under ideal conditions. Given a random variable x (in this case, log 
amplitude) with probability density / (x), mean , and variance o x 2 , we can select a random sample of 
size n and form the statistics X (sample mean) and S 2 (sample variance), where* 


_ 1 

X^~ 

n 



and 


( 4 ) 


S 2 



( 5 ) 


Both of these quantities are function of random variables and therefore are random variables also. 

With E being the expectation operator, 

E (X) = fi x (6) , • 


£(S2> = (^-‘) ^ ’ 

m 4 - m 2 2 2 (/n 4 -2m 2 2 ) m 4 -3 m 2 2 

E | [S 2 - E{S 2 )] 2 | = variance ( S 2 ) = — — + 

*The statistics X and S 2 are maximum likelihood estimators for n x and a x 2 . The amount of bias in S 2 becomes insignificant for 
for large samples. 


( 7 ) 

( 8 ) 
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where n is the number of independent samples, and m t is the i th central moment. Each data run 
consisted of 10,000 samples from the scintillation signal; however, because of the finite bandwidth of the 
scintillation, the number of independent samples is much less than 10,000. The number of independent 
samples is (Reference 5) 

n = 2 B eq T (9) 

where B eq is the equivalent white noise bandwidth of the signal, and T is the length of the time record; 
^eq is approximated by (Reference 5) 

n = ** 

eq 2 r°° R x (r) d t ( 10 > 

*/ — OO 

where R x (r) is the autocorrelation function for the scintillation at time-lag r. 

The authors have evaluated the autocorrelation function experimentally for both wavelengths and 
thereby determined the equivalent bandwidth. Typically, the bandwidth of the He-Ne signal is several 
times that of the C0 2 ; therefore, the dispersion of the He-Ne sample variance is much less than that of 
the C0 2 . In a strict sense, the equivalent bandwidth referred to in Equations 9 and 10 should be the 
bandwidth of the log irradiance signal and not the bandwidth of the irradiance fluctuations. However, 
the data reduction system did not permit an evaluation of the spectrum of log irradiance; therefore, the 
equivalent bandwidth of the irradiance signal was used. It is easy to show that the bandwidth of the 
irradiance signal is not seriously affected by the logarithmic operation under weak scintillation 
conditions (which is the situation at 10.6 microns, for most paths). Consider a system with input 
x = v 0 + v(t) and output y = A log x, 


where 

v 0 = constant (dc level), 

v(f) = random time function with mean value equal to zero, 

and 

A = constant (gain of the log-amplifier). 


Power spectral density (psd) analysis of the input x is performed on v(0 only, since v 0 merely results in a 
delta function at the origin in the frequency domain. The output of the system is 


y = A log [ Vq + v (01 = A log (v 0 ) 


1 + 


v(t) 


A log (v 0 ) + A log 


1 + 


v(0_ 

V n 


Therefore, psd analysis of the output needs to be performed only on the second term, since A log [v 0 ] is 
a constant. However, 
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log (1 + x) « X 


for x << 1 ; 


therefore, 


log 




for v(t) << v 0 (i.e., weak scintillation). 


Thus, psd analysis of y is equivalent to psd analysis of A v(t)/v 0 ; this is the same as psd analysis of the 
input x since the constant factors A, v 0 are removed by dividing through by the mean-square value. 

If n is known, and typical values for the log-amplitude central moments are used, numerical 
manipulation shows that, for a 1 -standard deviation scatter in the sample variances, the parameter i p can 
be expected to fluctuate by approximately ± 1.5. Therefore, the observed experimental scatter in i p is of 
the same order of magnitude as that which can be expected statistically from the finite data samples. 

The normalized power spectral density for both wavelengths for a particular data run is shown in 
Figure 5. These functions were evaluated by means of a Fourier transformation of their respective 
autocorrelation functions. Both curves are normalized by their intensity variances. Thus, if each 
wavelength has equal amounts of scintillation, this plot would indicate spectral distribution of the 
scintillation. The normalization scheme used requires the height of the C0 2 plot to be increased by a 
factor of 16. Therefore, it should be clear that, at any frequency, there is far more scintillation in the 

visible spectrum that at 10 microns. The 
resolution bandwidth is 5Hz, and the wind 
velocity perpendicular to the propagation path is 
approximately 3 mph. The power supply ripple 
on the 10-micron beam is obvious in this plot. 

Figure 6 presents the same data in the time 
domain in terms of the normalized 
autocorrelation function. Again, the internal 
modulation of the C0 2 laser appears as peaks at 
8- and 16-millisecond intervals. If one fits an 
exponential function of the form e~ aT to each of 
these curves, the correlation time constants (the 
alphas) have a ratio of 4.4, which is close to the 
ratio of the square root of the wavelengths. An 
alternate interpretation is that the time required 
for the autocorrelation function to decrease to 
the 1 /e point is proportional to the square root 
of wavelength. This experimental relationship 
has been noted in much of our data but has not 
been analyzed extensively. 

A typical portion of the analog scintillation 
signal is shown in Figure 7. The computer 



FREQUENCY (Hz) 


Figure 5— Dual wavelength normalized power spectral 
density (resolution bandwidth = 5Hz; each curve 
normalized to the same area; log-amplitude variances for 
these data are a 2 (Q 632) = 0.200 and ol (1Q 6) = 
0.00693). 
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r (milliseconds) 


Figure 6— Normalized auto correlation function for dual 
wavelength scintillation; fitting the data with the 
exponential e -or gives (C0 2 ) = 71.4, cdHe-Ne) = 312. 


evaluated the variances as 0.200 and 0.00668, 
which gives a \p value of 30. The chart recorder 
had a 1 10-Hz, 3-db bandwidth, which tends to 
smooth some of the high-frequency He-Ne 
scintillation (very little C0 2 scintillation is in 
the 1 10-Hz region). The important parameter is 
the ratio of the maximum and mininum signal 
excursions from the average value. For the 
He-Ne scintillation, excursions of several octaves 
are common. In both of these traces, the 
background levels were not significantly 
different from the zero level. 

Figure 8 is typical cumulative probability 
plot of log amplitude at 10.6 microns on 
Gaussian probability paper. The abscissa is F( &), 
where 


FU) 
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Figure 7— Typical analog scintillation signal; 
log-amplitude variances are 0.200 for He-Ne, and 
0.00668 for C0 2 . 



f(u ) being the probability density function for Figure 8— Cumulative probability of log amplitude at 

log amplitude. The ordinate is log amplitude, 10.6 microns; data plotted on a Gaussian abscissa scale, 

which is defined in Equation 1. The abscissa 
scale is constructed so that a normal random 

variable would plot as a a straight line. Thus it is clear that the statistics of scintillation at 10.6 microns 
satisfy this test for log-normality. 


7 



CONCLUSIONS 


In summary, the following conclusions can be made. 

1 Xhe average value for the log-amplitude variance ratio at 10.6 and 0.632 microns is very close 
to the predicted value of 26.7. 

2. The C0 2 scintillation has a bandwidth that is several times less than for the He-Ne, and a 
correlation time that is approximately 4.5 times greater than for the He-Ne. 

3. The probability distribution of scintillation at 10.6 microns is log-normal just as in the visible 
spectrum. 


Goddard Space Flight Center 

National Aeronautics and Space Administration 
Greenbelt, Maryland, November 14, 1968 
039-02-01-10-51 
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